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We give a short overview of the effective action approach in quan- 
tum field theory and quantum gravity and describe various meth- 
ods for calculation of the asymptotic expansion of the heat kernel for 
second-order elliptic partial differential operators acting on sections 
of vector bundles over a compact Riemannian manifold. We consider 
both Laplace type operators and non-Laplace type operators on mani- 
folds without boundary as well as Laplace type operators on manifolds 
with boundary with oblique and non-smooth boundary conditions. 
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1 Effective Action in Gauge Field Theories 
and Quantum Gravity 

In this lecture we briefly describe the standard formal construction of the 
generating functional and the effective action in gauge theories following 
the covariant spacetime approach to quantum field theory developed mainly 
by DeWitt [|l|. The basic object of any physical theory is the spacetime 
M, which is assumed to be a m-dimensional manifold with the topological 
structure of a cylinder 

M = /x S, 

where / is an open interval of the real line (or the whole real line) and S is 
some (m — l)-dimensional manifold. The spacetime manifold is here assumed 
to be globally hyperbolic and equipped with a (pseudo)-Riemannian metric 
g of signature ( — I — ■ +); thus, a foliation of spacetime exists into spacelike 
sections identical to S. Usually one also assumes the existence of a spin 
structure on M. A point x = {x^) in the spacetime is described locally by 
the time x^ and the space coordinates (x^, . . . , x"^~^). We label the spacetime 
coordinates by Greek indices, which run from to m — 1, and sum over 
repeated indices. 

Let us consider a vector bundle V over the spacetime M each fiber of 
which is isomorphic to a vector space, on which the spin group Spin(l, m — 1), 
i.e. the covering group of Lorentz group, acts. The vector bundle V can also 
have an additional structure on which a gauge group acts. The sections of the 
vector bundle V are called fields. The tensor fields describe the particles with 
integer spin (bosons) while the spin-tensor fields describe particles with half- 
odd spin (fermions). Although the whole scheme can be developed for the 
superfields (a combination of boson and fermion fields), we restrict ourselves 
in the present lecture to boson fields. A field if is represented locally by a 
set of functions 

if = (<^^(x)), 

where A = 1, . . . , dim V. Capital Latin indices will be used to label the local 
components of the fields. We will also use the condensed DeWitt notation, 
where the discrete index A and the spacetime point x are combined in one 
lower case Latin index 

i = {A,x). 



Ivan Avramidi: Heat Kernel Approach in Quantum Field Theory 2 

Then the components of a field (/? are 

As usual, we will also assume that a summation over repeated lower case 
Latin indices, i.e. a combined summation-integration, is performed, viz. 

J ■ Lf = Jiif^ = I dYo\{x) JA{x)ip (x), 
J M 

where dvo\{x) = dx\g\^^'^, \g\ = detg^^, is the natural Riemannian volume 
element defined by some background metric g. We will often omit the volume 
element when this does not cause any misunderstanding. 

In quantum field theory (QFT) the vector bundle V is called the con- 
figuration space. One assumes that the configuration space is an infinite- 
dimensional manifold Ai. The fields (f^ are the coordinates on this manifold, 
the variational derivative 6/6ip is a tangent vector, a small disturbance 6ip is a 
one-form and so on. If S{(p) is a scalar field on the configuration space, then 
its variational derivative SS/Sip is a one-form on Ai with the components 
that we denote by 

By using the functional differentiation one can define formally the concept 
of tangent space, the tangent vectors. Lie derivative, one-forms, metric, con- 
nection, geodesies and so on. 

The dynamics of quantum field theory is determined by an action func- 
tional S{ip), which is a differentiable real- valued scalar field on the config- 
uration space. The dynamical field configurations are defined as the field 
configurations satisfying the stationary action principle, i.e. they must sat- 
isfy the dynamical equations of motion 

6ip 

with given boundary conditions. The set of all dynamical field configurations, 
i.e. those that satisfy the dynamical equations of motion, A^o? is a subspace 
of the configuration space called the dynamical subspace. 

Quantum field theory is basically a theory of small disturbances on the 
dynamical subspace. Most of the problems of standard QFT deal with scat- 
tering processes, which are described by the transition amplitudes between 
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some well defined initial and final states in the remote past and the remote 
future. The collection of all these amplitudes is called the scattering matrix, 
or shortly S'-matrix. 

Let us single out in the space-time two causally connected in- and out- 
regions, that lie in the past and in the future respectively relative to the region 
Q, which is of interest from the dynamical standpoint. Let |in) and |out) be 
some initial and final states of the quantum field system in these regions. Let 
us consider the transition amplitude (out | in) and ask the question: how does 
this amplitude change under a variation of the interaction with a compact 
support in the region Q. The answer to this question gives the Schwinger 
variational principle which states that 

t 
6 (out I in) = — (out I 6S I in) , 
n 

where 6S is the corresponding change of the action. This principle gives 
a very powerful tool to study the transition amplitudes. The Schwinger 
variational principle can be called the quantization postulate, because all the 
information about quantum fields can be derived from it. 

Let us change the external conditions by adding a linear interaction with 
some external classical sources J in the dynamical region Q, i.e. 

6S = J ■ ip. 

The amplitude (out|in) becomes a functional of the sources that we denote 
by Z{J). By using the Schwinger variational principle one can obtain the 
chronological mean values in terms of the derivatives of the functional Z{J) 

(^^. ^n\ ^ (out|T(y-"---y-^)|in) 
^^ ^ / (out I in) 

i J 5Ji„ ■ ■ ■ 6Ji^ 

where T denotes the operator of chronological ordering that orders the (non- 
commuting) operators in order of their time variables from right to left. In 
other words the functional Z{J) is the generating functional for chronological 
amplitudes. Let us now define another functional W{J) by 
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We have obviously 




(^-. 


■■v)= -] e '^77 TT'' 


in particular, 






{0') = 0\ 
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5Ji^ ■■ -SJi 
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The functional is called the background or the mean field, ^ is called the 
full propagator, and ^**-*", are called the full connected Green functions or 
the correlation functions. Thus, whilst Z{J) is the generating functional for 
chronological amplitudes the functional W{J) is the generating functional for 
the connected Green functions. The Green functions satisfy the boundary 
conditions which are determined by the states |in) and |out). 

The mean field itself is a functional of the sources, (j) = 0(J), the derivative 
of the mean field being the full propagator 



■]t 



In the non-gauge theories the full propagator Q'^\ which plays the role of 

the (infinite-dimensional) Jacobian, is non-degenerate. Therefore, one can 

change variables and consider as independent variable and J(0) (as well as 

all other functionals) as the functional of 0. 

There are many different ways to show that there is a functional r(0) 

such that 

/6S{0)\_6T{<P) 

\ 60 / 5<p - 

This functional is defined by 

(outlin) = exp < — (T + J 

h 
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or by the functional Legendre transform 

r(0) = w{jm - j(0) ■ j^w{jm- 

This is the most important object in quantum field theory. It contains all 
the information about quantized fields. First of all, the first variation of F 
gives the effective equations for the background fields 



These equations replace the classical equations of motion and describe the 
effective dynamics of the background field with regard to all quantum correc- 
tions. That is why F is called the effective action. Furthermore, the second 
derivative of F(</)) determines the full propagator 

5^V 



The higher derivatives determine the so-called full vertex functions 

F 



,«i---«fc' 



which are also called strongly connected, or one-particle irreducible functions. 
In other words, F(0) is the generating functional for the full vertex functions. 
The full vertex functions together with the full propagator determine the full 
connected Green functions and, therefore, all chronological amplitudes and, 
hence, the S'-matrix. Thus, the entire quantum field theory is summed up in 
the functional structure of the effective action. 

One can obtain a very useful formal representation for the effective ac- 
tion in terms of functional integrals (called also path integrals, or Feynman 
integrals). A functional integral is an integral over the (infinite-dimensional) 
configuration space M.. Although a rigorous mathematical definition for the 
functional integrals is absent, they can be used in perturbation theory of QFT 
as an effective tool, especially in gauge theories, for manipulating the whole 
series of perturbation theory. The point is that in perturbation theory one 
encounters only the functional integrals of Gaussian type, which can be well 
defined effectively in terms of the classical propagators and vertex functions. 
The Gaussian integrals do not depend much on the dimension and, therefore. 
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many formulas from the finite-dimensional case, like the Fourier transform, 
integration by parts, delta-function, change of variables etc. are valid in the 
infinite-dimensional case as well. One has to note that the functional integrals 
are formally divergent — if one tries to evaluate the integrals, one encoun- 
ters meaningless divergent expressions. This difficulty can be overcome in 
the framework of the renormalization theory in so-called renormalizable field 
theories, but we will not discuss this problem in the present lectures. 

Integrating the Schwinger variational principle one can obtain the follow- 
ing functional integral: 

(out|in) = / Dip exp < - [S{ip) + J ■ ip]> . 
M 

Correspondigly, for the effective action one obtains a functional equation 
exp|^r(0)| = |l)v^exp|l 5M-^-(v^-0) }. 

M 

The only way to get numbers from this formal expression is to take advan- 
tage of the semi-classical approximation within a formal expansion in powers 
of Planck constant h: 

oo 

r = 5 + 5^n'^-ir(,). 

Substituting this expansion in the functional equation for the effective action, 
shifting the integration variable in the functional integral 

(p = (j) + vhh, 

expanding the action S{ip) in functional Taylor series in quantum fields h, 
expanding both sides of the equation in powers on h and equating the coef- 
ficients of equal powers of h, one gets the recurrence relations that uniquely 
define all coefficients ^{k)- All functional integrals appearing in this expansion 
have the form 

Vhexp (~-h-Ahj h'^ ■ ■ ■ h'" , 

M 

where A is a partial differential operator defined by the second variation of 
the action 

A-- — 
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These integrals are Gaussian and can be calculated in terms of the functional 
determinant of the operator A and the bare propagator G = A^^, i.e. the 
Green function of the operator A with Feynman boundary conditions, and 
the local classical vertex functions S^i^,„i^. 

In particular, the one-loop effective action is determined by the functional 
determinant of the operator A 

r(i) = --logDet A. 

Strictly speaking, the Gaussian integrals are well defined for elliptic dif- 
ferential operators in terms of the functional determinants and their Green 
functions. Although the Gaussian integrals of quantum field theory are de- 
termined by hyperbolic differential operators with Feynman boundary condi- 
tions they can be well defined by means of the analytic continuation from the 
Euclidean sector of the theory where the operators become elliptic. This is 
done by so-called Wick rotation — one replaces the real time coordinate by a 
purely imaginary one x^ -^ ir and singles out the imaginary factor also from 
the action S -^ iS and the effective action F -^ iV. Then the metric of the 
spacetime manifold becomes positive definite and the classical action in all 
'nice' field theories becomes a positive-definite functional. Then the fast os- 
cillating Gaussian functional integrals become exponentially decreasing and 
can be given a rigorous mathematical meaning. 

Let us try to apply the formalism described above to a gauge field theory. 
A characteristic feature of a gauge field theory is the fact that the dynamical 
equations 

are not independent — there are certain identities, called Nother identities, 
between them. This means that there are some nowhere vanishing vector 
fields 

R - R* J- 

on the configuration space M. that annihilate the action 

R,Q,iS = 0, 
and, hence, define invariance fiows on M.. The transformations of the fields 

5^ip' = R\C 
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are called the invariance transformations and R^ are called the generators of 
invariance transformations. The infinitesimal parameters of these transfor- 
mations C, are sections of another vector bundle (usually the tangent bundle 
TG of a compact Lie group G) that are respresented locally by a set of 
functions 

r) = in^)), 

a = l,...,dimG, over spacetime with compact support. To distinguish 
between the components of the gauge fields and the components of the gauge 
parameters we introduce lower case Latin indices from the beginning of the 
alphabet; the Greek indices from the beginning of the alphabet are used as 
condensed labels 

a = (a, x) 

that include the spacetime point. We assume that the vector fields R^ are 
linearly independent and complete, which means that they form a complete 
basis in the tangent space of the invariant subspace of configuration space. 
The vector fields Rq, form the gauge algebra. We restrict ourselves to the sim- 
plest case when the gauge algebra is the Lie algebra of an infinite-dimensional 
gauge Lie group Q, which is the case in Yang-Mills theory and gravity. Then 
the flow vectors Rq decompose the configuration space into the invariants 
subspaces of Ai (called the orbits) consisting of the points connected by the 
gauge transformations. The space of orbits is then M./Q. The linear inde- 
pendence of the vectors R^ at each point implies that each orbit is a copy 
of the group manifold. One can show that the vector fields R^ are tangent 
to the dynamical subspace A^o, which means that the orbits do not intersect 
M.Q and the invariance flow maps the dynamical subspace M.q into itself. 
Since all field configurations connected by a gauge transformation, i.e., the 
points on an orbit, are physically equivalent, the physical dynamical variables 
are the classes of gauge equivalent field configurations, i.e., the orbits. The 
physical configuration space is, hence, the space of orbits M./Q. In other 
words the physical observables must be the invariants of the gauge group. 

To quantize a gauge theory by means of the functional integral, we con- 
sider the in- and out- regions, define some |in) and |out) states in these 
regions and study the amplitude (out | in). Since all field configurations along 
an orbit are physically equivalent we have to integrate over the orbit space 
Ai/Q. To deal with such situations one has to choose a representative field in 
each orbit. This can be done by choosing special coordinates {I^{ip),x°'{^)) 
on the configuration space Ai, where I^ label the orbits and x" the points 
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in the orbit. Computing the Jacobian of the field transformation and intro- 
ducing a delta functional 6{x — C) ^^ can fix the coordinates on the orbits 
and obtain the measure on the orbit space M./Q 

where 

is a non-degenerate operator. Thus we obtain a functional integral for the 
transition amplitude 

(out|in) = J^V^DetFi^)5ixiv) " C) exp {^^m} • 

Now one can go further and integrate this equation over parameters ( with a 
Gaussian measure determined by a nondegenerate matrix 7 = (70/3), which 
most naturally can be choosen as the metric on the orbit (gauge group met- 
ric). As a result we get 



(out|in) = / r'v'Det^/^7DetF(v9)exp 



M 



s{v) + ^x{v) -ixiv) 



The functional equation for the effective action takes the form 



exp||r(0)l = fvipDet^/^-fBetFiip) 



M 



X exp 

h 



's{v) + ^x(v^) ■ ixiv) - ^ ■ (v^ - <^)^ 



This equation can be used to construct the semi-classical perturbation theory 
in powers of the Planck constant (loop expansion), which gives the effective 
action in terms of the bare propagators and the vertex functions. The new 
feature is though that the bare propagator and the vertex functions are de- 
termined by the action 

Scsiv) = S{^) + -x{<^) ■ ixiv) 

h h 

+ -logDetF(^) + — logDet7 
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In particular, one finds the one-loop effective action 

1 1 1 

Tfi) = -— logDetA + -logDetF + — logDet7, 

where 

2 Heat Kernel Asymptotic Expansion 

As we have seen in the previous lecture the effective action in quantum field 
theory can be computed within the semi-classical perturbation theory — the 
one-loop effective action is determined by the functional determinants of 
second-order hyperbolic partial differential operators with Feynman bound- 
ary conditions and the higher-loop approximations are determined in terms 
of the Feynman propagators and the classical vertex functions. As we noted 
above these expressions are purely formal and need to be regularized and 
renormalized, which can be done in a consistent way in renormalizable field 
theories. One should stress, of course, that many physically interesting the- 
ories (including Einstein's general relativity) are perturbatively non-renor- 
malizable. Since we only need Feynman propagators we can do the Wick 
rotation and consider instead of hyperbolic operators the elliptic ones. The 
Green functions of elliptic operators and their functional determinants can 
be expressed in terms of the heat kernel. That is why we concentrate in the 
subsequent lectures on the calculation of the heat kernel. 

The gauge invariance (or covariance) in quantum gauge field theory and 
quantum gravity is of fundamental importance. That is why, a calculational 
scheme that is manifestly covariant is an inestimable advantage. A manifestly 
covariant calculus is such that every step is expressed in terms of geometric 
objects; it does not have some intermediate "non- covariant" steps that lead 
to an "invariant" result. Below we describe a manifestly covariant method 
for calculation of the heat kernel following mainly our papers [0, |^, ^, ^]. 

Let (M, g) be a smooth compact Riemannian manifold of dimension m 
without boundary, equipped with a positive definite Riemannian metric g. 
Let l^ be a vector bundle over M, V* be its dual, and End {V) = V ®V* he 
the corresponding bundle of endomorphisms. Given any vector bundle V", we 
denote by C^iV) its space of smooth sections. We assume that the vector 
bundle V is equipped with a Hermitian metric. This naturally identifies the 
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dual vector bundle V* with V, and defines a natural L^ inner product and 
the L^-trace Tr 1^2 using the invariant Riemannian measure on the manifold 
M. The completion of C°°{V) in this norm defines the Hilbert space L'^{V) 
of square integrable sections. We denote by TM and T*M the tangent and 
cotangent bundles of M. Let a connection, V^ : C°°{V) -> C°°{T*M V), 
on the vector bundle V be given, which we assume to be compatible with 
the Hermitian metric on the vector bundle V. The connection is given its 
unique natural extension to bundles in the tensor algebra over V and V*. 
In fact, using the Levi-Civita connection V^*" of the metric g together with 
V^, we naturally obtain connections on all bundles in the tensor algebra over 
V, V*, TM and T*M; the resulting connection will usually be denoted just 
by V. It is usually clear which bundle's connection is being referred to, from 
the nature of the section being acted upon. Let V* be the formal adjoint to V 
defined using the Riemannian metric and the Hermitian structure on V and 
let Q e C°°(End {V)) be a smooth Hermitian section of the endomorphism 
bundle End(\/). 

A Laplace type operator F : C°°{V) -^ C°°(y) is a partial differential 
operator of the form 

F = V*V + Q = -g^'V^V, + Q. (1) 

It is obviously symmetric, i.e. 

iFip,^P) = {ip,F^|J), 

elliptic, and can be made essentially self-adjoint, i.e. its closure is self-adjoint, 
which implies that it has a unique self-adjoint extension. We will not be very 
careful about this and will simply say that F is elliptic and self-adjoint. It 
is well known that: 

i) the operator F has a discrete real spectrum, {Xn}'^=i, bounded from 
below: 

A < Ai < A2 < ■ ■ ■ < A„ < ■ • • 

with some real constant A, 

ii) all eigenspaces of the operator F are finite-dimensional, and iii) the 
eigenvectors, {(Pn}'^=i, of the operator F, are smooth sections of the 
vector bundle V that form a complete orthonormal basis in L'^[V). 
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For t > the operators 

U{t) = exp(-tF) 

form a semi-group of bounded operators on L'^{V), so called heat semi-group. 
The kernel of this operator is defined by 

U{t\x,x') = J2e~'^"Vn{x)®V*^{x'), 
n 

where each eigenvalue is counted with multiplicities. It is a section of the 
external tensor product of vector bundles V Kl y* over M x M, which can 
also be regarded as an endomorphism from the fiber of V over x' to the fiber 
of V over x. This kernel satisfies the heat equation 

{dt + F) U{t) = (2) 

with the initial condition 

U(Q-^\x,x) = 5{x,x) (3) 

and is called the heat kernel. 

Moreover, the heat semigroup U{t) is a trace-class operator with a well 
defined L^-trace 

Tr ^2 exp(-tF) = f tiv^'^^it) . (4) 

M 

Hereafter try denotes the fiber trace and the label 'diag' means the diagonal 
value of a two-point quantity, e.g. 

f/dias(t|a;) = U{t\x,x') 

The trace of the heat kernel is obviously a spectral invariant of the operator F. 
It determines other spectral functions by integral transforms. Of particular 
importance is the so-called zeta function, which enables one to define, in 
particular, the regularized functional determinant of an elliptic operator. 

In these lectures we will study the heat kernel only locally, i.e. in the 
neighbourhood of the diagonal of M x M, when the points x and x' are close 
to each other. We fix a point x' of the manifold and consider a small geodesic 
ball with a radius smaller than the injectivity radius of the manifold, so that 
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each point x of the ball can be connected by a unique geodesic with the 
point x'. Such geodesic ball can be covered by a single coordinate patch with 
normal local coordinates centered at x'. Let a{x, x') be the geodetic interval, 
defined as one half the square of the length of the geodesic connecting the 
points X and x', i.e. 

(7{x,x') = -[dist (Xjx')]^. 

The first derivatives of this function with respect to x and x' define tangent 
vector fields to the geodesic at the points x and x' 



u''{x,x') = g'^'d^a, 



u" 



'{x,x') = g^''''dl,a, (5) 



and the determinant of the mixed second derivatives defines the so-called 
Van Vleck-Morette determinant 



1 , , „ , 1 



A{x,x') = \g{x)\-^\g{x')\"det{-d,d',,a) . (6) 

Let, finally, V{x, x') denote the parallel transport operator of sections of the 
vector bundle V along the geodesic from the point x' to the point x. It is 
an endomorphism from the fiber of V over x' to the fiber of V over x (or 
a section of the external tensor product V Kl V* over M x M). Near the 
diagonal of M x M all these two-point functions are smooth single-valued 
functions of the coordinates of the points x and x' . We should point out from 
the beginning that we will construct all two-point geometric quantities (in 
particular, the coefficients of the asymptotic expansion of the heat kernel as 
t — > 0) in form of covariant Taylor series. The Taylor series do not necessarily 
converge in smooth case; they do, however, converge in the analytic case in 
a sufficiently small neighborhood of the diagonal. 
Further, one can easily prove that the function 



[/o(t) = (47rt)— /W/2 exp (-^) P 



satisfies the initial condition (|^). Moreover, locally it also satisfies the heat 
equation (H) in the free case, when the Riemannian curvature Riem of the 
manifold, the curvature TZ of the connection V^, and the endomorphism Q 
vanish: 

Riem = 7^ = Q = 0. 
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Therefore, t/o(t) is the exact heat kernel for the Laplacian in flat Euclidean 
space with a flat trivial bundle connection. This function gives a good frame- 
work for the approximate solution in the general case. Namely, by factorizing 
out the free factor we get an ansatz 

Uit) = (4vrt)-™/2Ai/2exp (-^) pf)(t) . (7) 

The function Q{t\x,x'), called the transport function, is a section of the 
endomorphism bundle End {V) over the point x'. It satisfies the transport 
equation 



(dt + jD + l\ nit) = 0, 



with the initial condition 

n{o\x,x') = 1, 

where I is the identity endomorphism of the bundle V (we will often omit it) 
and D and L are differential operators defined by 

D = u^V„ (8) 

L = V-^A-^'^FA''^V. (9) 

Now, let us fix a sufficiently large negative parameter A, viz. A < A, so 
that {F — AI) is a positive operator. Since 

exp[-t(F - AI)] = e*^exp[-tF], 

the transport function for the operator {F — AI) is e^^VL{t). Clearly, for suffi- 
ciently large negative A, A << 0, the function e*'^f2(t) with all its derivatives 
decreases faster than any power of t as t -^ oo. Let us consider a slightly 
modified version of the Mellin transform of the function e*'^i7(t) 

1 /"°° 

^^^^^ = f{^J, dtr'~'e''n{t). (10) 

Note that for fixed A this is a Mellin transform of e*'*'r2(t) and for a fixed 
q this is a Laplace transform of the function t~'^~^Q{t). The integral (|10D 



converges for Reg < 0. By integrating by parts and analytical continuation 
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one can prove that the function bg{X) is an entire function of q. The values 
of the function bq{X) at the integer positive points q = k are given by 

n=0 ^ ^ 

where 

ak = i-dt)''nit) . (12) 

t=o 

These coefficients a^ = ak{x,x') are called Hadamard-Minakshisundaram- 
DeWitt-Seeley (HMDS) coefficients. 

By inverting the Mellin transform we obtain a new ansatz for the trans- 
port function and, hence, for the heat kernel 

c+ioo 



^(t) = ^ / dqe-'H'^T{-q)b,{X), 



where c < and Re A < A. Clearly, since the left-hand side of this equation 
does not depend on A, neither does the right hand side. Thus, A serves as an 
auxiliary parameter that regularizes the behavior at t ^ oo. 

Substituting this ansatz into the transport equation we get a functional- 
differential equation for the function bg 



(l + -D\bgi\) = iL-\I)bg_^i\) 



with the initial condition 

bo{X) = I. 

Note that for integer q = k and A = this becomes a differential recursion 
system for the coefficients Uk 

ao = I, (13) 

1 + -Djak = Lak^i. (14) 

It is interesting to note that there is an asymptotic expansion of bq{\) as 

A -^ — oo 

^n!r(g-n + l) 
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that coincides with (|TT]) for integer q. 

By computing the inverse Melhn transform we obtain the asymptotic 
expansion of the transport function as t — > in terms of the coefficients ak 

A:=0 

Using our ansatz (|^ we also find the trace of the heat kernel in form of 
an inverse Mellin transform 

c+joo 

Tr^2exp(-tF) = (47rt)— /V*^-i^ / rfgt'^r(-g)i?,(A), 

ITXI J 
c—ioo 

where 

B,{X) = ltTvbf^^{X). (16) 

M 

Noting that Bg is an entire function of q, this gives the standard asymptotic 
expansion as t ^ 

oo 

Tr L^ exp(-tF) ~ Yl ^^'""^^'^a. , (17) 



fc=0 



where 
and 



A2k+i = (18) 

^. = (47r)-"/2H)!y"trvaf^ (19) 

M 

This is the famous Minakshisundaram-Pleijel asymptotic expansion, which is 
called Schwinger-DeWitt expansion in the physics literature. This expansion 
is of great importance in differential geometry, spectral geometry, quantum 
field theory and other areas of mathematical physics, such as theory of Huy- 
gens' principle, heat kernel proofs of the index theorems, Korteveg-De Vries 
hierarchy, Brownian motion etc.. 

The (off-diagonal) HMDS coefficients a^ are determined by the recursion 
system (|T^ . The formal solution of this recursion system is 

Uk = D;' L D^\ L---D^^ LI, (20) 
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where 

Dk = l + ^D. (21) 

To give a precise meaning to this formal operator solution we need to 
define the inverse operator D^^. This can be done in terms of the covariant 
Taylor series. We will need the following notions from the theory of symmet- 
ric tensors. Let S*^ be the bundle of symmetric tensors of type (m, n). First 
of all, we define the exterior symmetric tensor product 

\/ • C" v C* V cn+i 

V . D„ X Dj- -^ '->m+j 

of symmetric tensors by 

This naturally leads to the following definition of the exterior symmetric 
power of a symmetric tensor 

\/k . qn qnk 

V . 0„j ^mk 

^ (23) 

k 

Next, we define the inner product 

, . qn qi qi 

We also define the exterior symmetric covariant derivative 

VS cm qm 

by 

(V^A)J;;t^^=V(„,Aj;:f:^^). (25) 

These definitions are naturally extended to End (V) -valued symmetric ten- 
sors, i.e. to the sections of the bundle S^^ End {V). 
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Let us consider the space of smooth two-point functions in a small neigh- 
borhood of the diagonal x = x' that we will denote by |/). Let us define a 
special set of such functions |n), labeled by a non-negative integer n, by 

|0) = 1, (26) 

i-lY 
\n) = L^^v^m', (27) 

n'. 

where u' is the tangent vector field to the geodesic connecting the points x 
and x' at the point x' defined by (|]). It is easy to show that these functions 
satisfy the equation 

D\n)=n\n) (28) 

and, hence, are the eigenfunctions of the operator D with positive integer 
eigenvalues. 

Let (n| denote the dual functions defined by 

W) = {V'rf , (29) 

x=x' 

so that 

(n|m) = (5mnl(n) , (30) 

where I(^n) is the identity endomorphism on the space of symmetric ra-tensors. 
Using this notation the covariant Taylor series for an analytic function |/) 
can be written in the form 

\f) = Y,\n)^{n\f), (31) 

n=0 

and, therefore, the functions \n) form a complete orthonormal basis in the 
subspace of analytic functions. 

The complete set of eigenfunctions \n) can be employed to present the 
action of the operator L on a function |/) in the form 

L\f)= E \m)^{m\L\n)^{n\f), (32) 

m,n>0 

where {m\L\n) are the 'matrix elements' of the operator L that are just 
End (\^)-valued symmetric tensors, i.e. sections of the vector bundle ^^(M)® 
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End (V). Now it should be clear that the inverse of the operator D^ in (2T) 
can be defined by 

oo , 

^^•'l/) = EiT^H^H/)- (33) 

Using such representations for the operators D^^ and L we obtain a covariant 
Taylor series for the coefficients ak 

oo 

a'k = '^\n)-k{n\ak) (34) 

n=0 

where 

{n\ak) = Y. ri-r^M^I^I^^-i) 

'k{nk-i\L\nk-2) -k- ---k (?2i|L|0) , (35) 

with Uk = n. 

Thus, we have reduced the problem of computation of the HMDS-coefficients 
Ofc to the evaluation of the matrix elements ijn\L\n) of the operator L. For a 
differential operator L of second order, the matrix elements {m\L\n) vanish 



for n > 171 + 2. Therefore, the summation over rii in ( p5D is limited from 
above: rii > 0, and Ui < nj+i + 2, for z = 1, 2, . . . , fc — 1, and, hence, the sum 
([35| ) always contains only a finite number of terms. 

We will not present here explicit formulas, (they have been computed 
explicitly for arbitrary m, n in [^), but note that all these quantities are 
expressed polynomially in terms of three sorts of geometric data: 

i) symmetric tensors of type (2, n), i.e. sections of the bundle 5*^ obtained 
by symmetric derivatives 

Kin) = (V^)"-2Riem 

of the symmetrized Riemann tensor Riem taken as a section of the 
bundle S'|, 

n) sections 
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of the vector bundle S'^®End {V) obtained by symmetrized derivatives 
of the curvature TZ of the connection V^ taken as a section of the bundle 

Si ^ End (V), 

iii) End (K)-valued symmetric forms, i.e. sections of the vector bundle 
S^ End {V), constructed from the symmetrized covariant derivatives 

Qin) = (V^)"Q 

of the endomorphism Q. 

3 Approximation Schemes for Calculation of 
the Heat Kernel 

In this lecture we are going to investigate the general structure of the heat 
kernel coefficients Ak. We will follow mainly our papers |0, |^, ||, ^, ^ (see 
also our review papers [|ll| |l^, |l^). Our analysis will be again purely local. 
Since locally one can always expand the metric, the connection and the endo- 
morphism Q in the covariant Taylor series, they are completely characterized 
by their Taylor coefficients, i.e. the covariant derivatives of the curvatures, 
more precisely by the objects -ft'(n), 7^(n) and (5(n) defined in the previous 
lecture. We introduce the following notation for all of them 

and call these objects covariant jets; n will be called the order of a jet 3fJ(n). It 
is worth noting that the jets are defined by symmetrized covariant derivatives. 
This makes them well defined as ordering of the covariant derivatives becomes 
not important — it is only the number of derivatives that plays a role. 

The coefficients A^ are integrals of local invariants tryo^'*^^ which are 
polynomial in the jets. The first two coefficients have the well known form 

at' = I, (36) 

af^ = Q-U, (37) 

o 

where R is the scalar curvature. For k >2 one can classify the terms in A2k 

according to the number of the jets and their order 

k 
A2k = ^^2fc,(n) , 
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where A2kJn) can be presented symbolically in the form 



^2A:,(2) = / try 2_^^{0)^{ 



M 



/2k-6 
try 2_^ 2_^3fi(o)3ff(j)3fJ 
™- — n 



2k-6 

(2fc-6-i); 



M ^=0 






M ^=0 



/trvE^J 



^2A:,(A:) — / try ^ ^^(Q)- 



M 

More precisely, all quadratic terms can be reduced to a finite number of 
invariant structures, viz. ||^ 

^2.,(2) = (4vr) (-1) 2(2A;_3)! 



Ik 

M 



+fi'^R% a '-'R^. + ffR n '^'rI, 



where n = g^'^'VfjV^, and fj: are some numerical coefficients. These numer- 
ical coefficients can be computed by the technique developed in the previous 
section. The result reads 

Jk ~ -^' 



fk 



(2) ^ 1 

2{2k- 1) 
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A3) 

h 


K — L 


2{2k-iy 


/(4) 
Jk 


1 


2(4F - 1) 


f(5) _ 
Jk 


k'^-k-l 


4(4A;2 - 1) 



One should note that the same results were obtained by a completely different 



method in 1113 



Let us consider the situation when the curvatures are small but rapidly 
varying, i.e. the derivatives of the curvatures are more important than the 
powers of them. Then the leading derivative terms in the heat kernel are the 
largest ones and, therefore, the trace of the heat kernel has the form 

Tri2exp(-tF) = t-^/^^o + ^'"'"^^2 

where H{t) is some complicated nonlocal functional that has the following 
asymptotic expansion as t — ;> 



H{t)^^t''"^A 



2k,{2)- 



k=2 



Using the results for A2kX2) one can easily construct such a functional H just 
by a formal summation of leading derivatives 



H{t) = (47r) 



-m/2 



M 



+27^%Vc.-7^'^-^^)V/37^^^ - 2Q7(3)(-t n)i? 
+i?-^^{4) (_i n)i?^, + i?7^'^ i-t n)R\ , (38) 

where Rap is the Ricci tensor, and 7^*^(2;) are some entire functions defined 
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byi 

1 



where 



7«(^)= I dU^HOewl-^-^z 



f'\0 = 1, 



K3) 



b 

Therefore, H{t) can be regarded as generating functional for quadratic terms 
^2fc,(2) (leading derivative terms) in all coefficients A2k- It plays a very im- 
portant role in investigating the nonlocal structure of the effective action in 
quantum field theory in so-called high-energy approximation ||^. 

Let us consider now the opposite case, when the curvatures are strong 
but slowly varying, i.e. the powers of the curvatures are more important 
than the derivatives of them. Since the derivatives are naturally identified 
with the momentum (or energy), a situation when the derivatives are small is 
called the low-energy approximation in physics literature. The investigation 
of the low-energy effective action is of great importance in quantum gravity 
and gauge theories because it describes the dynamics of the vacuum state of 
the theory. The main terms in this approximation are the terms without any 
covariant derivatives of the curvatures, i.e. the lowest order jets. We will 
consider mostly the zeroth order of this approximation which corresponds 
simply to covariantly constant background curvatures 

VRiem = 0, V7^ = 0, VQ = 0. 

The trace of the heat kernel has then the form 

Tr i2 exp(-tF) = f-'^/^Qit) + 0{V^) , 
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where G(t) is a functional that has the following asymptotic expansion as 

oo 
A:=0 

where A2k,{k) are the terms without covariant derivatives (highest order terms 
in the jets) in the coefficients A2k and 0(V3fJ) denotes the terms with at least 
one derivative that vanish in the covariantly constant case. The terms A2k,k 
are just polynomials in the curvatures and the endomorphism Q. Therefore, 
the functional 0(t) is a generating functional for all heat kernel coefficients 
Ak for a covariantly constant background, in particular, for all symmetric 
spaces. 

There is a very elegant indirect way to construct the heat kernel without 
solving the heat equation but using only the commutation relations of differ- 
ential operators ^ ^, y, |TU|. The main idea is in a generalization of the usual 



Fourier transform to the case of operators and consists in the following. Let 
us consider for a moment a trivial case, where the curvatures vanish but the 
potential term does not: 

Riem = 0, 7^ = 0, VQ = 0. 

In this case the operators of covariant derivatives obviously commute and 
form an Abelian Lie algebra, i.e. 

It is easy to show that the heat semigroup operator can be presented in the 
form 

exp(-tF) = (47rt)-"/2gxp(-tQ) 

X f dk\g\'/'exp(-^-^ + k-vY 

where {k, k) = k^g^yk^ and A; ■ V = k^^V ^. Here, of course, it is assumed that 
the covariant derivatives also commute with the metric, i.e. 

[V,^] = 0. 
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Acting with this operator on the Dirac distribution and using the obvious 
relation 

exp{k-V)6{x,x') =S{k), (39) 

x=x' 

one integrates easily over k and obtains the diagonal of the heat kernel and 
the trace 

Tri2exp(-tF) = (47rt)-™/2 Ar y exp(-tQ). 

M 

Let us consider now a more complicated case when there is a nontrivial 
covariantly constant curvature 7^ 7^ in fiat space: 

Riem = 0, V7^ = 0, VQ = 0. 

In this case the covariant derivatives form a nilpotent Lie algebra 

[v^,v.] = 7^^., 

For this algebra one can prove a theorem expressing the heat semigroup 
operator in terms of an average over the corresponding Lie group 

exp(-tF) = (4^^)"™/'e^P(-^Q)det^(!i(™)(-^) 



X 
R 



/*l9l"=exp 



— {k,tncoth{tn)k) 



exp {k ■ V) 



Here functions of the curvature 7^ are understood as functions of sections of 
the bundle End (TM) ® End(\^), and the determinant det End (tm) is taken 
with respect to End {TM) indices, End {V) indices being intact. 

It is not difficult to show that in this case the equation (|39|) is still valid, 
so that the integral over k^ becomes trivial and we obtain immediately the 
trace of the heat kernel 

Tr^.exp(-tF) = (4vrt)--/2 |tr y exp(-tg)det ^/„^,(^,,) (— ^) . 

M 



Ivan Avramidi: Heat Kernel Approach in Quantum Field Theory 26 



Expanding this in a power series in t one can find all covariantly constant 
terms A^ ^ in all heat kernel coefficients A^. 

Let us now generalize the algebraic approach to the case of curved man- 
ifolds with covariantly constant Riemann curvature and trivial connection 

VRiem = 0, 7^ = 0, VQ = 0. 

First of all, we give some definitions. The condition VRiem = de- 
fines the geometry of locally symmetric spaces. A Riemannian locally sym- 
metric space which is simply connected and complete is a globally sym- 
metric space (or, simply, symmetric space). A symmetric space is said to 
be of compact, noncompact or Euclidean type if all sectional curvatures 
K{u, v) = Rap'^su^^v^u^v^ are positive, negative or zero. A direct product 
of symmetric spaces of compact and noncompact types is called semisimple 
symmetric space. A generic complete simply connected Riemannian sym- 
metric space is a direct product of a flat space and a semisimple symmetric 
space. 

It should be noted that our analysis is purely local. We are looking for a 
universal (in the category of locally symmetric spaces) local generating func- 
tion of the curvature invariants, that reproduces adequately the asymptotic 
expansion of the trace of the heat kernel. This function should give all the 
terms without covariant derivatives of the curvature A2k,{k) in the asymptotic 
expansion of the heat kernel, i.e. in other words all heat kernel coefficients 
A2k for any locally symmetric space. It turns out to be much easier to obtain 
a universal generating function of t whose Taylor coefficients reproduce the 
heat kernel coefficients A2k than to compute them directly. 

It is obvious that no flat subspace contributes to the coefficients A2k- 
Therefore, to find this universal structure it is sufficient to consider only 
semisimple symmetric spaces. Moreover, since the coefficients Ak are poly- 
nomial in the curvatures, one can restrict oneself only to symmetric spaces 
of compact type. Using the factorization property of the heat kernel and 
the duality between the compact and the noncompact symmetric spaces one 
can obtain then the results for the general case by analytical continuation. 
That is why we consider only the case of compact symmetric spaces when 
the sectional curvatures and the metric are positive definite. 

First of all, we choose a basis for the tangent bundle TM that is parallel 
(covariantly constant) along the geodesies. The frame components of the 
curvature tensor of a symmetric space are, obviously, constant and can be 
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presented in the form 



R. 



abed 



(^ikE\bE cd, 



where E^at, i = 1, ■ ■ ■ ,p, is some set of antisymmetric m x m matrices, with 
p being a constant satisfying p < m{m — l)/2, and f3ik is some symmetric 
p X p nondegenerate matrix. The traceless matrices Di = (-D'^*) defined by 

D°'ib = —f3ikE cbg^"" = —E>"'bi 

are known to be the generators of the holonomy algebra Ti 

[D,,Dj,]=F\kD,, 

where F-'-. are the structure constants. 



ik 



In symmetric spaces a much richer algebraic structure exists. Indeed, let 
us define the quantities C^^bc ~ '^^ 
by 



-C^CBi A = ^1 ■ ■ ■ iD, where D = m + p, 



-C/ ab) ^ lb 



jja 



lb} 

/^a /^i /^a 

•-^ be — '-^ ka — '-^ ik 

and the matrices Ca = {C^ac) = {Ca,Ci): 

D'a 



= 0, 



kl 



F' 



kh 



Ca 



a 



EK 



^^ la 

F^ 



ik 



One can show that they satisfy the Jacobi identities p, |T0 

[Ca, Cb\ = C abCc 



and, hence, define a Lie algebra Q of dimension D with the structure con- 
stants C'^BC, the matrices Ca being generators of the adjoint representation. 
In symmetric spaces one can find explicitly the generators of the infinites- 
imal isometrics, i.e. the Killing vector fields ^a, and show that they form a 
Lie algebra of isometrics that is (in case of semisimple symmetric space) 
isomorphic to the Lie algebra Q, viz. 



[a,e 



B 



C ABic ■ 
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Moreover, introducing a symmetric nondegenerate D x D matrix 



lAB 



gab 

/3ik 



that plays the role of the metric on the algebra Q, one can express the operator 
F in semisimple symmetric spaces in terms of the generators of isometrics 



where '-y^^ = {'Jab)'^- 

Using this representation one can prove that the heat semigroup operator 
can be presented in terms of an average over the group of isometrics G [^, |T^ 



exp(-tF) 



(47rt) 



-D/2 



exp 



-t{Q- -Rg 






ikhl^'^det\'Lj''''f^^f) 



X exp 



Ad(e?) 



--{k,k) + k-^ 



where I7I = det7AB, {k, k) = k'^'jABk^, k ■ C = k^CA, k ■ ^ = k'^^A, and 
Rg is the scalar curvature of the group of isometrics G 



R, 



G 



ABr^C 



■iD 



~t7 ^ AD^ BC- 



Acting with this operator on the Dirac distribution 6{x, x') one can, in 
principle, evaluate the off-diagonal heat kernel exp{—tF)6{x, x'), i.e. for non- 
coinciding points X ^ x' (see fl^l)- To calculate the trace of the heat kernel, 
it is sufficient to compute only the coincidence limit x = x'. Splitting the 
integration variables k^ = (g°, cj*) and solving the equations of characteristics 
one can obtain the action of the isometrics on the Dirac distribution f0|, |TU 

'smh{u-D/2) 



exp {k ■ ^) 6{x, x') 
where u ■ D = u'^Dj. 



dpf ~^ 

^^^ End (TM) 



uj-D/2 



5{q). 
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Using this result one can easily integrate over q to get the heat kernel 
diagonal. After changing the integration variables u -^ \/tuj it takes the 



form P, 10 



^diag^^^ = (47rt)-"^/2 



exp 



— 1 1 Q R Rh 

' ^ 8 6 






■{u,uj) 



1/2 ( smh{^tuj-F/2) \ 
^^^^Am)\ ^tu.F/2 J 



xdet 



F/2 
_i/2 f smh{Viuj-D/2) \ 

End (TM)^ V^^.D/2 J 



(40) 



where /? = det/5jj, (d^jCj) = u^PijU^, u ■ F = u^Fi, Fi = (F^ik) are the 
generators of the holonomy algebra 7i in adjoint representation and 



1 



^H — — tP -T ji-r km 



is the scalar curvature of the holonomy group. 

The remaining integration over u in (^Op can be done in a rather formal 



way ||T^, |T3l . Let a* and a^ be some operators acting on a Hilbert space that 
form the following p-dimensional Lie algebra 

[a\al] =51 

[a\a''] = [alal] =0. 
Let |0) be the 'vacuum vector' in the Hilbert space, i.e. 

(0|0) = 1, 

a'\0) = 0, 
(0|< = 0. 



Ivan Avramidi: Heat Kernel Approach in Quantum Field Theory 



30 



Then the heat kernel (^D|) can be presented in an algebraic form without any 
integration, i.e. 



f/diag^^^ 



(47rt) 



-m/2 



exp 



-t ( Q — o -^ ~ c -^^ 



x(0 



det 



1/2 



f sinh {Via ■ F/2) 



^^(^) V Via . F/2 



^*^^^ End (TM) 



sinh {Vta ■ D/2) 



Via ■ D/2 



X exp((a*,/3 ^a*)) 



where a ■ F = a F^. and a ■ D = a D^. This formal solution should be 
understood as a power series in the operators a'^ and a^; it determines a well 
defined asymptotic expansion as t ^ 0. 

By expanding these formulas in an asymptotic power series as t — *> one 
obtains all HMDS-coefficients a^!^^ for any locally symmetric space. Thereby 
one finds all covariantly constant terms A2k,{k) in all heat kernel coefficients. 



4 Heat-kernel Asymptotics for Non-Laplace 
Type Operators 

In this lecture, we study a general class of second-order non-Laplace type 
elliptic partial differential operators, acting on sections of a vector bundle 
V over a Riemannian manifold M without boundary following our papers 
15], |16| . In general, the study of spin-tensor quantum gauge fields in a gen- 



eral gauge necessarily leads to non-Laplace type operators acting on sections 
of general spin-tensor bundles described in the first lecture. It is precisely 
these operators that are of prime interest in the present lecture. The study 
of non-Laplace operators is quite new, and the available methods are still 
underdeveloped in comparison with the Laplace type theory. The only ex- 
ception to this is the case of anti-symmetric forms, which is pretty simple 
and, therefore, is well understood now [|l^, [l^, |l^, ^, ^ . 

So, we will restrict our attention to operators acting on tensor-spinor 
bundles. These bundles may be characterized as those appearing as direct 
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summands of iterated tensor products of the cotangent and spinor bundles, 
i.e. V = TM (g) ■ ■ ■ (g) TM (g) T*M ® ■■■ T*M ® 5, with S being the spinor 
bundle. Alternatively, they may be described abstractly as bundles asso- 
ciated to representations of the spin group Spin(m). These are extremely 
interesting and important bundles, as they describe the fields in Euclidean 
quantum field theory. The connection on the tensor-spinor bundles is built in 
a canonical way from the Levi-Civita connection. The generators are deter- 
mined by the representation of Spin(m) which induces the bundle V; they are 
tensor-spinors constructed purely from Kronecker symbols, together with the 
fundamental tensor-spinor if spin structure is involved. More general bundles 
appearing in field theory are actually tensor products of these with auxiliary 
bundles, usually carrying another (gauge) group structure. 

Let Q be a smooth Hermitian section of the bundle End {V), i.e. Q* = Q, 
and a be a parallel symmetric Hermitian End (V^)-valued tensor, more pre- 
cisely, a smooth section of the vector bundle TM (gTM® End (V^) satisfying 
the following conditions 



a"" = a""" 



Va = 0. (41) 

The operator of our primary interest in this lecture has the form 

F = V*aV + Q 

= -a'^^V^V. + g. (42) 

Non-Laplace type operators appear naturally in the context of Stein- Weiss 
operators [^. Let 

T*M (S)V = Wi®---®Wn 

be the decomposition of the bundle T*M ^V in irreducible components Wj 
and 

Pr^ : T*M 0V ^Wi 

be the corresponding projections. Stein- Weiss operators Gi : C°°{V) —>■ 
C°°{Wi) are first order partial differential operators (called the gradients) 
defined by 

Gi = Pr,;V. 
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Then the operator F : C^{V) -^ C^{V) defined by 

n 

with some constants q, is a second-order non-Laplace type operator of the 



form (^2]) with 

n 

a = J]Q(Pr,rPri. 

Now it is obvious that the structure of the coefficient a depends solely on the 
representation of the spin group with which the bundle V is associated. 

We will restrict ourselves to a special class when the coefficient a is built 
in a universal, polynomial way, using tensor product and contraction from 
the metric g and its inverse g*, together with (if applicable) the volume form 
e and/or the fundamental tensor-spinor 7. Such a tensor-endomorphism a is 
obviously parallel. Here we do not assume that a^^'^ has the form g^'^ly oi" 
g^^'^B with some automorphism B G Aut iV) . We do not set any conditions 
on the endomorphism Q, except that it should be Hermitian. 

In the following we will denote the leading symbol of the operator F, 
o"l(F;x,^), with ^ a cotangent vector, just by A{x,^). For the non-Laplace 
type operator F in (^) it has the form 

We require that the leading symbol should be positive definite, i.e. A{x, ^) is 
a Hermitian and positive definite endomorphism for any (x, ^), with ^ 7^ 0. 
In particular, F is elliptic. Positive definiteness implies that the roots of the 
characteristic poljTiomial 

x(a;,e,A):=detv'[A(a;,0-A] 

are positive functions on M. 

A very important point is that the structure of the spectrum, i.e. the 
number s of eigenvalues, Ai, . . ., A^, and their multiplicities di , . . . ,ds are 
constant on M. Moreover, one can show that tryy4"'(a;, ^), and, therefore, 
the characteristic polynomial x(a;,^. A) depends on (x, ^) only trough |^p = 
g^^{x)^^^u- As a result, the dependence of the eigenvalues Ai(x,^) on (x, ^) 
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is only through |^p as welL Since A{x, ^) is 2-homogeneous in ^, the Aj must 
be also: 

for some positive real numbers fii , . . . , Hs, which are independent of the point 
(x,^) G T*M, and, in fact, independent of the specific Riemannian manifold 
{M,g). 

Let IIj be the orthogonal projection onto the Aj -eigenspace. The IIj 
satisfy the conditions 



{t^k) 



n? = 


= n,. 


n^Efc = 


= 0, 


s 

>;n. = 

i=l 


= ly, 


tryHi = 


= d,. 



In contrast to the eigenvalues, the projections depend on the direction ^/|^| of 
^, rather than on the magnitude |^|. In other words, they are 0-homogeneous 
in ^. Furthermore, they are even polynomials in ^/|^|: 

P 1 

n=0 1*^1 

where p is some positive integer. Here the Ili{2n) are some End (V)-valued 
trace- free symmetric 2n-tensors that do not depend on ^. Clearly, the leading 
symbol can be written in terms of eigenvalues and projections 

s 

^(a^,0 = ierE/^^n,(x,o. (43) 

There is also a converse formula for the projections in terms of powers of 
the leading symbol. Note that the highest degree of projections, p, is also a 
constant that depends only on the representation to which V is associated; 
both s and p can be computed explicitly in representation-theoretic terms 



16|. 



The non-Laplace type operator F with a positive definite leading symbol 
is an elliptic self-adjoint operator of second order. Therefore, there is a well 



Ivan Avramidi: Heat Kernel Approach in Quantum Field Theory 34 

defined heat kernel U{t\x,x'). Moreover, there is a well defined heat kernel 
diagonal f/*^'^^(t) and the trace of the heat kernel Tr j;^2 exp(— tF) that have 
the asymptotic expansion of the standard form (|1^. Since the global heat 
kernel coefficients Ak are determined by the integrals of the fiber trace of the 
local ones a^^,"^^, it is sufficient to compute the local heat kernel coefficients, 
more precisely, their fiber traces, tr^a^'^^. By invariance theory, these coeffi- 
cients are linear combinations of the local invariants built from the geometric 
objects (curvatures, the potential Q, and their derivatives) with universal nu- 
merical constants. It is these universal constants that we want to compute. 
Therefore, this can be done at any fixed point of the manifold. 

Let us stress here that our purpose is not to provide a rigorous con- 
struction of the heat kernel with estimates; for this we rely on the standard 
references 0. Rather, given that the existence of heat kernel asymptotic 
expansion is known, our aim is to compute its coefficients. 

Our analysis will be again purely local. We fix a point x' in the manifold 
M and consider a small geodesic ball with the radius smaller than the injec- 
tivity radius of the manifold. Then any point in this ball can be connected 
with the fixed point x' by a unique geodesic. Further, we represent the heat 
kernel in the form 

U{t\x,y) = A^/^{x,x')V{x,x')U{t\x,y;x') 

xV-\y,x')A'/\y,x'), (44) 

where A is the Van Vleck-Morette determinant and V is the parallel transport 
operator defined in lecture 2. Then the modified heat kernel W is a section 
of the bundle End (V) at x' but is scalar at x and y. It satisfies the modified 
heat equation 

(dt + L)U{t) = (45) 

where L = V^^ A^^^"^ F A^^'^V is the operator defined by (^, with the initial 
condition 

V({0'^\x, y; x') = A~^{x, x')5{x, y) . 

Here and everywhere below, as usual, the differential operators act on the 
first space argument of the heat kernel (recall that x' is being fixed). 

We shall employ the standard scaling device for the heat kernel lA{t\x, y; x') 
when a; — >■ x', 1/ — > x', and t — > 0. We introduce a small expansion parameter 
£, choose the normal coordinates at x' , and scale the coordinates according 
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to 



X ^- x + e(x — x) 



y -^ 


x' + e{y - x') , 




t -^ 


e\ 


(46) 


that this also means that the derivatives scale according to 




dt 


^ ^ 
^ i^^*' 




9, 


^ \o.. 





Note also that in normal coordinates 

^{x,x') = \g{x)\'^/\ 

so that 

A"^(x, x)5{x, y) = 6{x - y), 

and 

V{x,x') = I 

(however, VP ^ 0!). 

Next, we expand the operator L and the heat kernel W in a formal asymp- 
totic series in e 

oo 
n=0 

and 



W(t)~5^e"-'"f/„(t). 

ra=0 

The zeroth order heat kernel is determined by the equation 

{dt + Fo)Uo{t) = (47) 

with the initial condition 

Uo{0+\x,y;x') = 6{x-y). (48) 
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The higher order approximations are determined by the following differential 
recurrence relations 

fe-i 



idt + Fo)Uk{t) = -Y,Fk-nUn{t) 



n=0 

with the initial conditions 

[/fc(0+|x,y;x') = 0. 
By construction the coefficients Un are homogeneous functions, i.e. 

[/„(t|a;,y;x') = t^"-'")/^ x f/„ ("f |x' + ^-^,x' + ^^; x'V (49) 

Therefore, on the diagonal one obtains the asymptotic expansion 

oo 
n=0 

where 

Uf''^{t\x)=Un{t\x,X;x). 

Comparing this with the standard heat kernel asymptotic expansion ([T7|) we 
see that the diagonal odd-order coefficients vanish 

and the even-order ones give the heat kernel coefficients 

M 

Using the form of the operator L in the leading order 

Fo = -a^-[x')d^d, (51) 

we easily find the leading order heat kernel by Fourier transform 

f/o(tk,i/;x') = j ^e^«-(-2/)exp[-tA(x',0], 

R™ 
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where ^ ■ (x — y) = $,^{x'^ — y^). Here and everywhere below all integrals over 
^ will be over R™. Writing the leading symbol in terms of the projections, 
we get 

The trace of the diagonal can now be easily computed 

s 

tr vUtm = 5^rf.(47rt/i,)''"/^ • (52) 

It gives the diagonal value of the lowest order heat kernel coefficient oq : 

+^ ^diag _ V^ "i 
tryCg - 7 ^ 

and, therefore, 

Ao=(4vr)-/^X:4,vol(M). 

1=1 f^i 

These formula points out a new feature of non-Laplace type operators; one 
which complicates life somewhat. Whereas the dimension dependence of the 
heat coefficients of Laplace type operators is isolated in the overall factor 
(47r)^™/^, the dimension dependence for non-Laplace type operators is more 
complicated. 

The calculation of higher-order coefficients is a challenging task. We will 
indicate how the coefficient Ai is computed. By the invariance theory we 
have, 



_ m/2 



A2= ftTv{HQ + IpR) , 



M 

where /9 is a universal constant and H is some endomorphism. Both j3 and 
H depend only on the leading symbol of the operator F. 

To compute these quantities we need, first of all, the Taylor expansion of 
the metric and connection in normal coordinates 
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+0[{x-x')% 

Q{x) = Q + 0[{x-x')]. (53) 

Here and below all coefficients are computed at the fixed point x'; we do 
not indicate this explicitly. Similarly, the Taylor expansion of the tensor- 
endomorphism a^^^{x) is determined by the equation V^a"^ = 0, which gives 

a^"'{x) = a^"' + -a^^^R''^aX(3ix''-x"'){x^-x'^) 

+0[{x-x'f]. 
Using these formulas we obtain 
Fi = 0, 
F2 = X''\p{x''-x"'){x''-x"')d^d, 

where 

X^^aP = -o"^ R ia\X\l3), 
2 1 

and [A, -B]+ = AB + BA denotes the anticommutator. 

From the recurrence relations (and the initial conditions) we find 

and 



U2{t) = -JdTUo{t-T)F2Uo{r) 
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By using the Fourier representation for Uq we obtain for the diagonal 

t 



where 

F2 = X^\pdpli,i, - Y\d'^i^ + Q . 

This gives finally the heat kernel coefficient A2 

1 



A2 = -(47r)-™/2 11^ /"drtrye-(i-^)^(^'«)F2e-^^("'«). 



M 

It is not very difficult to obtain from here the endoniorphism if; it is given 
by 

The calculation of the coefficient /? is a much more complicated problem. 
After a long calculation one obtains a complicated expression in terms of the 
constant /ij (the eigenvalues of the leading symbol) and the leading symbol 
of the operator F (see, [p!5[] ). 

An interesting feature of the non-Laplace type operators is the semi- 
classical polarization. This means that, unlike Laplace type operators, the 
asymptotic expansion of the heat kernel for non-Laplace type operators has 
the form 



cr(x,a;') 



2t/ii 



VLi{t\x,x') 



U{t\x,x') = ^(47rt/i,)-"^/2^i/2(^^^/)g^p 

where each transport function Vti satisfies a different transport equation. This 
also implies that the differential recursion system for the coefficients of the 
asymptotic expansion of the transport functions is much more complicated. 

5 Heat-Kernel Asymptotics of Oblique 
Boundary- Value Problem 

In this lecture we study the heat kernel asymptotics for a Laplace type par- 
tial differential operator acting on sections of a vector bundle over a compact 
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Riemannian manifold with boundary. In this case one has to impose some 
boundary conditions in order to make a (formally self-adjoint) differential 
operator self-adjoint (at least symmetric) and elliptic. There are many ad- 
missible boundary conditions that guarantee the self-adjointness and ellip- 
ticity of the problem. The simplest boundary conditions are the classical 
Dirichlet and the Neumann ones. There exist also slight modifications of the 
Neumann boundary conditions (called Robin boundary conditions in physical 
literature) when the normal derivative of the field at the boundary is pro- 
portional to the value of the field at the boundary [^, |2^ . In an even more 



general scheme, called mixed boundary conditions, the Dirichlet and Robin 
boundary conditions are mixed by using some projectors. In this lecture 
we study a more general setup, called the oblique boundary-value problem, 
which includes both normal and tangential derivatives of the fields at the 
boundary [E31 EH|. We will follow mainly our papers |20, 



Let (M, g) be a smooth compact Riemannian manifold of dimension m 
with smooth boundary dM with a positive-definite Riemannian metric g on 
M and induced metric g on dM. In this lecture the Greek indices range from 
1 through m and lower case Latin indices range from 2 through m. Let 

{ei} = {e2,...,em}, 

be the local frame for the tangent bundle TdM and x = (x*) = (x^, . . . , x™"), 
be the local coordinates on dM. Let r be the normal geodesic distance to 
dM, and 

be the inward pointing unit normal to dM. Let K be a (smooth) vector 
bundle over the manifold M, V be a connection on V and V be the induced 
connection on the boundary. Further let Q be a smooth endomorphism of V 
and F be a Laplace type operator 

F = -g'"'V,V, + Q. 

Let W = V\qj^ he the restriction of the vector bundle V to the boundary 
dM. We define the boundary data map ^ : C°°{V) ^ C°°{W ® W) by 

Let n be an orthogonal Hermitian projector acting on W and F be an anti- 
Hermitian End (ly)-valued vector field F on the boundary that satisfies the 
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conditions 






-r, 
rn = 0. 



Let, further, S he a smooth Herniitian endoniorphisni of the bundle W or- 
thogonal to n, i.e. 



S* 



US 



s 



^n = o. 



We use these to define a first-order self-adjoint tangential differential operator 
A: C'^{W)^C'^{W) by 



A = (I - n) <^{rv, + v,r) + ^1 (I - n) . 

We study the oblique boundary-value problem for the Laplace type op- 
erator F. The oblique boundary conditions now read 

where B : C°°{W © W) -^ C°°(W © W) is the boundary operator defined by 



B 



n 
A i-n 



(54) 



This is equivalent to the following boundary conditions 



(I - U)Vn^ 



dhl 



Yiif 



+ Ay? 



dM 



dM 



0. 



(55) 



The boundary operator B (|S3) incorporates all standard types of boundary 
conditions. Indeed, by choosing 11 = I and A = one gets the Dirichlet 
boundary conditions, by choosing 11 = 0, A = I one gets the Neumann 
boundary conditions. More generally, the choice T = corresponds to the 
mixed boundary conditions. 
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Integration by parts shows that the Laplace-type operator F endowed 
with obhque boundary conditions is symmetric, i.e. 

However, it is not necessarily elliptic. To be elliptic the boundary-value 
problem {F, B) has to satisfy two conditions. First of all, the leading symbol 
of the operator F should be elliptic in the interior of M. Second, the so-called 
strong ellipticity condition should be satisfied. This question was studied in 
p6| where it has been shown that the oblique boundary-value problem for 
a Laplace type operator is strongly elliptic with respect to the cone C \ R+ 
if and only if for any nonvanishing cotangent vector C, on the boundary the 
endomorphism If^lI — iT-C is positive-definite, i.e. |C|I — "JT-^ > 0. A sufficient 
condition for strong ellipticity is: 

|Cri + (r-c)'>o. 

The heat kernel U{t\x,y) is now defined by the heat equation 

(9i + F)f/(t) = 0, 

the initial condition 

f/(0+|x,?/) =5{x,y) , 

and the boundary conditions 

Bij[U{t\x,y)]=^. (56) 

Hereafter the boundary data map (as well as the boundary operator) acts on 
the first argument of the heat kernel. 

The heat kernel of a smooth boundary- value problem on a manifold with 
boundary has the following asymptotic expansion as t — >■ 0"*" 

oo 

Tr L2 exp{-tF) ~ ^ t^'^-™)/^^^ . (57) 

fc=0 

In the case of manifolds without boundary only even order terms were present 
(see p!7|) . Now, in contrast, all A^ are non-zero. They have the following 
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general form: 



4 _ / ^(0) , / „(l) 
^2fc — / a,i, + / 0- 



2k - I "'2k ^ j "2fc' 
M dM 



^2k+i - / a^k+i 

dM 



Hereafter the integration over the boundary is defined with the help of the 
usual Riemannian volume element dvol g on dM with the help of the induced 
metric g. 

Here a^ are the (local) interior heat-kernel coefficients and a^ are the 
boundary ones. The local interior coefficients a^ are determined by the 
same local invariants as in the manifolds without boundary, i.e. by the 
HMDS coefficients trv-a^'^^, and therefore, do not depend on the boundary 
conditions. The boundary coefficients a^ are far more complicated because 
in addition to the geometry of the manifold M they depend essentially on the 
geometry of the boundary dM and on the boundary conditions. For Laplace- 
type operators they are known for the usual boundary conditions (Dirichlet, 
Neumann, or mixed version of them) up to Og |]2^ , 23|. For oblique boundary 



conditions including tangential derivatives some coefficients were recently 



computed in [^, ^, ^ |29[. In this lecture we will evaluate the coefficient 



Ai, following our recent work pB|] . 

Let us fix a positive number S > 0. We split the whole manifold in a 
disjoint union of two different parts: 

M = M'''^ U M^^^^ 

where M^^"^ is a narrow geodesic strip near the boundary dM of the width 
5 and Af^* is the interior of the manifold M (without the thin strip), i.e. 
M'^* = M\ M^""^. 

We will construct the parametrix on M by using different approximations 
in M^^'^ and M™*. Strictly speaking, to glue them together in a smooth way 
one should use 'smooth characteristic functions' of different domains (parti- 
tion of unity) and carry out all necessary estimates. What one has to control 
is the order of the remainder terms in the limit t — > and their dependence 
on 6. Since our task here is not to prove the form of the asymptotic expan- 
sion, which is known, but rather to compute explicitly the coefficients of the 
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asymptotic expansion, we will not worry about such subtle details. We will 
compute the asymptotic expansion as t — * in each domain and then take 
the limit 5 — »• 0. 

We will use different local coordinates in different domains. In M™* we 
can, for example, choose normal coordinates centered at a fixed point xq. 
In fact, this is not necessary — we can use a manifestly covariant technique 
described in lecture 2. In M^'^'^ we choose the local coordinates as follows. 
By using the geodesic flow we get the local frame 

{N,e,} 

for the tangent bundle TM and the local coordinates x = {r,x) on M^'^'^, 
which identifies M^""^ with dM x (0,5). 

The construction of the parametrix in the interior f/™*(t|x, y) goes along 
the same lines as for manifolds without boundary described in the previous 
lectures. The idea is always to separate the basic case (when the coefficients 
of the operator F are frozen at a fixed point Xq). In the case of manifolds 
without boundary the basic case is, in fact, zero-dimensional, i.e. algebraic. 
The interior parametrix is defined by the heat equation (|^), the initial condi- 
tion dl) and by an asymptotic condition at infinity (instead of the boundary 
conditions). This means that effectively one introduces a small expansion 
parameter e reflecting the fact that the points x and y are close to each other 
and the parameter t is small. This can be done by fixing a point Xq = x' 
in M'"*, choosing the normal coordinates at this point (with gfj,u{x') = 5^u)-i 



scaling like in ( P6| ) and expanding in a power series in e. This construction is 
standard and we do not repeat it here (see lectures 2 and 4). Locally, at any 
point in M™*, the resulting interior parametrix is given by the same formulas 
as the heat kernel for a manifold without boundary, i.e. by the same formulas 
as in the lecture 2. For a fixed finite 5 > the error of this approximation 
is exponentially small as t — > 0. Thus, the interior parametrix has the same 
asymptotic expansion with t — > as the heat kernel for a manifold without 
boundary. In other words, as t — > 

/ trvf7^i(t) ~ f;t^-/2(4vr)-/2H)! f tr yaf^^ , (58) 

where Uf^^^^ are the standard local HMDS coefficients for manifolds without 
boundary computed in lecture 2. By taking the limit 6 ^ oi this equation 
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we obtain 

4ti = 0, (59) 

4? = (47r)-™/2H)!trv^af^ (60) 

For an elliptic boundary-value problem the diagonal of the parametrix 
U^^{t) in M^'^'^ has exponentially small terms, i.e. of order ~ exp(— r^/t), 
as t ^ 0+ and < r < S. These terms behave like distributions near 
the boundary, and, therefore, the integrals over M*^"*^, more precisely, the 
integrals 

lim||rfr(...), 

dM 

do contribute to the asymptotic expansion with coefficients being the inte- 
grals over dM. It is this phenomenon that leads to the boundary terms in 
the heat kernel coefficients. Thus, such terms determine the local boundary 
contributions a^ to the global heat-kernel coefficients A^. 

The boundary parametrix U^'^'^{t\x,x') in M'^'^'^ is constructed as follows. 
Now we want to find the fundamental solution of the heat equation near 
diagonal, i.e. for x ^ x' and for small t -^ in the region M^^'^ close to the 
boundary, i.e. for small r and r', that satisfies the boundary conditions on 
dM and an asymptotic condition at infinity. We fix a point on the bound- 
ary xq G dM and choose normal coordinates on dM at this point (with 
gij{0,Xo) =Sij). 

To construct the boundary parametrix, we again scale the coordinates. 
But now we include the coordinates r and r' in the scaling 

X ^ xq + e{x — xq), 

X -^ Xq + b{x — Xo) 

r -^ er, 
r' -^ er', 

_2 



eH. (61) 
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The corresponding differential operators are scaled by 

d ^ -9, 

e 

e 
dt ^ -,dt. (62) 

Then, we expand the scaled operator F in a power series in e, i.e. 

oo 
?i=0 

where -F„ are second-order differential operators with homogeneous symbols. 
Next, we expand the scaled boundary operator (with an extra factor e at the 
operator A) 

oo 
n=0 

where -B(„) are first-order tangential operators with homogeneous symbols. 
At zeroth order we have 

Fo = -dl - d\ 

Ho 

-Do = I 

where Ho = liix') and 



^°^\Ao I-Ho 



92 = g^\x')d,dk, 

Ao = T^{x)dj. (63) 

Note that all leading-order operators Fq, Bq and Aq have constant coefficients 
and, therefore, are very easy to handle. 

The subsequent strategy is rather simple. We expand the scaled heat 
kernel in e 



00 

J2—m+n J jhnd 

n 1 

ra=0 



f/bnd _ Y^ e2-m+nfjt 
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and substitute into the scaled version of the heat equation and the boundary 
condition. Then, by equating the hke powers in e one gets an infinite set of 
recursive differential equations for f/„ 



(9, + Fo)f/r = -5^F„f/r„, 



n=l 



with the boundary conditions 



B,i,[Ut-^] = -Y.Bni'iUt^ 



n=l 



and the asymptotic condition at infinity 

lim Ul'"^{t\r,x]r\x') = 0, 



(64) 



In other words, we decompose the parametrix into the homogeneous parts 
with respect to {x — Xq), {x' — Xq), t, r' and t. By using this homogeneity 
we obtain finally the asymptotic expansion of the diagonal of the boundary 
parametrix 



^dtjw~E^^'~™^^'^^'(i 



A;=0 



r 



Vi' ' Vt 



Now we have to integrate the diagonal f/^iag o'^^^ M^^^, expand it in an 
asymptotic series as t ^ 0, and then take the limit 5^0. One should stress 
that the volume element rfvol {x) = y^\g\dx should also be scaled, i.e. 



oo ^ 



dvo\ {er, x) ~ rfvol (0, x) ■ >^ ^^-rrdk 



fc=0 



X) 



where 



9k[x) 



Qk 



dvo\ (r, x) 



r=0 



(ivol (0,a;) 

The coefficients Qk will contribute directly to the coefficients of the asymptotic 
expansion. 
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We have 







dr 



M^nd 



dM 



X 



dvol (r, x) 



tivU^'^U 1 



Vt'^^'Vi'' 



dvol (0, x) 

Since as 5 ^ the volume of M^^*^ vanishes, i.e. 

hmvol(M^°'^) = 0, 

the contribution of all regular terms will vanish in the limit S ^ 0. In 
contrary, the singular terms, which behave like distributions near dM, will 
give the dM contributions to the boundary heat kernel coefficients a^ . By 

<5 

changing the integration variable r = \/i^ the integral f dr{. . .) becomes 





S/Vi 



dit 



1/2/ 



and in the limit t ^ becomes the improper integral 



dit 



1/2/ 



plus an exponentially small remainder term. Then in the limit 5 -^ we 
obtain integrals over DM up to an exponentially small function that we are 
not interested in. More precisely, as the result we get the coefficients a\. in 
the form 



.(1) 



fc-i -^/^ 

X^l^.iim /rfer 



xtrvf/n-i(l|e,5:;e,£). 



(65) 
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In this lecture we will find the boundary parametrix of the heat equation 

-bi 




to leading order, i.e. Uq^'^. We fix a point x' G dM on the boundary and the 



normal coordinates at this point (with gik{x') = 
TOM and replace the manifold M by 

Mn = TdM X R , . 



5ik), take the tangent space 



By using the explicit form of the zeroth-order operators Fq, Bq and Aq we 
obtain the equation 



and the boundary conditions 



0, 



i-Ti,){dr + indAu^''%t\x,y) 



r(a;)=0 



r(a;)=0 



(66) 



(67) 



(68) 



where Ho = n(x'), Fq = V^{x'). As usual the differential operators always act 
on the first argument of a kernel. Moreover, for simplicity of notation, we will 
denote Ho and Fq just by FI and F-' and omit the dependence of all geometric 
objects on x' . To leading order this does not cause any misunderstanding. 
Furthermore, the heat kernel should be symmetric and vanish at infinity. 

By using the Laplace tarnsform in t and Fourier transform in [x — y) this 
equation reduces to an ordinary differential equation of second order in r on 
R+, which can be easily solved taking into account the boundary conditions 
at r = and r ^ oo. Omitting simple but lengthy calculations we obtain 



Urit\x,y) 



dC 



W + lOO 



R" 



(27r) 



m— 1 



d\ 

2Td 



xe 



-tX+i(-{x- 



y^G{X\C,r{x),r{y)), 



(69) 



where w is a negative constant and G is the leading-order resolvent kernel in 
momentum representation. It reads 



G{X\C,u,v) 



1 



^vW^ 
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xl exp l-\u - v\^/\C\^ - X\ 

I - 2n + 2^r ■ c (VICP-A - ^r ■ c 



-1 



X exp 



-(^ + ^)v¥F^A 



where Re ^KP - A > 0. 

By changing the integration variables, deforming the contour of integra- 
tion and computing certain Gaussian integrals, we obtain the heat kernel 
diagonal 



U^'^'itlr^x^r^x 



) = (47rt)-'"/2|i 



+ exp,-!^](I-2n) + c,^^ 



where 



$(z) 



dC 



R"i-i C 



duj 



X 



exp[-|cp-cj2 + 2itu2] r-c(cui + r-c)~^ 



(70) 



Here the contour of integration C comes from — oo + ie, encircles the point 
uj = i\(\ in the clockwise direction and goes to +oo + ie, with £ > a positive 
infinitesimal parameter; the contour C does not cross the interval Recu = 0, 
< Imuj < \(\, on the imaginary axis and is above all singularities of the 
resolvent G. 

Now by using eq. (|65|) we obtain the coefficient a\ 



,(1) 



oo 

(4vr)-™/2 l^^tr V. {e-«'(I - 2n) + $(0} 



and, finally, by computing the integral over ^ we get 

.1 



a? 



(4t) 



~{m~l)/2 ■ 



i-2n 
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+2 / ;^exp[-|Cr^-(r.C)T}. (71) 

We now consider two particular cases. First of all, if the matrices F* form 
an Abelian algebra, i.e. 

[r,F^] = o, 

then the integral ( [TTD is Gaussian and can be easily evaluated explicitly: 

a« = (47r)-("'-i)/2i |_i _ 2n + 2(1 + T^^/^} . 

Another very important case is when the operator A is a natural Dirac type 
operator when the matrices F-' form a Clifford like algebra 

pp + rip = 2g'^- — ^ — -F^, 
[m — 1) 

where F^ = g^jT^T^ . In this case one obtains 

a« = (47r)-(™-i)/2l|_j_2n 



[m — 1) 



-(m-l)/2 



Note that the integral (^) diverges when the strong ellipticity condition, 
ICpI + (F ■ C)^ > 0, is violated. This leads to singularities in the heat kernel 
coefficients. This is a general feature of the oblique boundary- value problem. 

6 Heat-Kernel Asymptotics for Non-Smooth 
Boundary Conditions 

The boundary-value problem studied in the previous lecture was in the 
smooth category. A more general (and much more complicated) setting, so 
called singular boundary- value problem, arises when either the symbol of the 
differential operator or the symbol of the boundary operator (or the bound- 
ary itself) are not smooth. In this lecture we study a singular boundary- value 
problem for a second order partial differential operator of Laplace type when 
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the operator itself has smooth coefficients but the boundary operator is not 
smooth. The case when the manifold as well as the boundary are smooth, but 
the boundary operator jumps from Dirichlet to Neumann on the boundary, 
is known in the literature as Zaremba problem. Zaremba problem belongs to 
a much wider class of singular boundary- value problems, i.e. manifolds with 
singularities (corners, edges, cones etc.). There is a large body of literature 
on this subject where the problem is studied from a very abstract function- 
analytical point of view (see |3y] and the references therein.) However, the 



study of heat kernel asymptotics of Zaremba type problems is quite new, and 
there are only some preliminary results in this area ||31], ^ . Moreover, com- 
pared to the smooth category the needed machinery is still underdeveloped. 
In this lecture we will closely follow our papers [ P5| , ^^. 

Let (M, g) be a smooth compact Riemannian manifold of dimension m 
with smooth boundary dM with a positive-definite Riemannian metric g on 
M and induced metric g on dM. In this lecture we will be dealing with 
submanifolds of Riemannian manifolds of codimension one and two. There- 
fore, we need to fix notation, ffist of all. With our notation, Greek indices, 
IJ,,^, . . ., label the local coordinates on M and range from 1 through m, lower 
case Latin indices from the middle of the alphabet, i,j,k,l,..., label the local 
coordinates on dM (codimension one manifold) and range from 2 through m, 
and lower case Latin indices from the beginning of the alphabet, a,b,c,d, . . ., 
label the local coordinates on a codimension two manifold Sq C dM that 
will be described later and range over 3, . . . , m. Further, we will denote by 
g the induced metric on the submanifolds (of the codimension one or two). 
We should stress from the beginning that we slightly abuse the notation by 
using the same symbols for all submanifolds (of codimension one and two). 
This should not cause any misunderstanding since it is always clear from the 
context what is meant. 

Let V^ be a vector bundle over the manifold M, V be a connection on 
V and V be the induced connection on the boundary. Further let Q be a 
smooth endomorphism of V and F be a Laplace type operator 

F = -g^'V^V, + Q. 

Let W = V^lg^-f be the restriction of the vector bundle V to the boundary 
dM and N be the inward pointing unit normal to the boundary. We use 
these to define the boundary data map ijj : C°°{V) -^ C°°{W © W) and 
the boundary operator B : C°°{W ® W) ^ C°°{W ® W); the boundary 
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conditions then are 

Btpiip) = 0. 

We always assume the manifold M itself and the coefficients of the op- 
erator F to be smooth in the interior of M. If, in addition, the boundary 
dM is smooth, and the boundary operator 5 is a differential operator with 
smooth coefficients, then {F, B) is a smooth local boundary- value problem. 
If the boundary dM consists of a finite number of disjoint connected parts, 

dM = U^^iS,, 

with each Sj being compact connected manifold without boundary, 

d^i =0 

and 

Ei n Ej = 

if i 7^ J, then one can impose different boundary conditions on different 
connected parts of the boundary Ej. This means that the full boundary 
operator decomposes 

fi = fii © • • ■ © 5„, 

with Bi being different boundary operators acting on different bundles. If 
the boundary operators are smooth (even if different), then such a boundary- 
value problem is still smooth. 

In this lecture we are interested in a different class of boundary conditions. 
Namely, we do not assume the boundary operator to be smooth. Instead, 
we will study the case when it has discontinuous coefficients. Such problems 
are often also called mixed boundary conditions; to avoid misunderstanding 
we will not use this terminology. We impose different boundary conditions 
on connected parts of the boundary, which makes the boundary-value prob- 
lem discontinuous. Roughly speaking, one has a decomposition of a smooth 
boundary in some parts where different types of the boundary conditions 
are imposed, i.e. say Dirichlet or Neumann. The boundary operator is then 
discontinuous at the intersection of these parts. We consider the simplest 
case when there are just two components. We assume that the boundary of 
the manifold DM is decomposed as the disjoint union 

aM = El U E2 U Eo , (72) 
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where Si and S2 are smooth compact sub manifolds of dimension {m — 1) 
(codimension 1 submanifolds) , with the same boundary 

So = 9Si = 9E2, 

that is a smooth compact submanifold of dimension {m — 2) (codimension 2 
submanifold) without boundary, i.e. 

9So = 0. 

Let us stress here that when viewed as sets both Si and S2 are considered 
to be disjoint open sets, i.e. 

SinS2 = 0. 

Let TTi and tt2 be the trivial projections of sections, ip, of a vector bundle 
W to Sj defined by 

if X G Sj and 

if a; ^ Sj. In other words tti maps smooth sections of the bundle W to 
their restriction to Si, extending them by zero on S2, and similarly for tt2. 
Let S G C°°{End{W)) be a smooth Hermitian endomorphism of the vector 
bundle W. 

We study the Zaremba type boundary- value problem for the Laplace type 
operator F. The Zaremba boundary conditions are 

where B : C°^{W(BW) — > C°^{W(BW) is a Zaremba type boundary operator 
defined by 

\^7r267r2 712 J 

In other words, we have Dirichlet boundary conditions on Si and Neumann 
(Robin) boundary conditions on S2: 

V =0, (74) 



{VN + S)ip = 0. (75) 
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In the following, for simplicity, we restrict ourselves to the case S = 0. The 
projectors tti and 7T2 as well as the boundary operator B are clearly non- 
smooth (discontinuous) on Eq. Note that the boundary conditions are set 
only on open subsets Si and E2; the boundary conditions do not say anything 
about the boundary data on Sq. 

By integrating by parts on dM, it is not difficult to check that the 
Zaremba type boundary- value problem [F, B) for a Laplace-type operator 
with the boundary operator B ( ffSD is symmetric. One can also show that it 
is elliptic with respect to C \ R+. 

The heat kernel is defined by the equation 

{dt + F)U{t\x,x') = 0, (76) 

with the initial condition 

U{0+\x,x')=6{x,x'), (77) 

and the boundary condition 

Bij[U{t\x,x')]=0. (78) 

Since coefficients of the boundary operator B are discontinuous on Sq, 
Zaremba type boundary- value problem is essentially singular. For such prob- 
lems the asymptotic expansion of the trace of the heat kernel has additional 
non-trivial logarithmic terms pO[, i.e. 



00 



TTL2exp{-tFB) ~ 5^t('^-")/25fc + logt^t'=/2i/fc. (79) 

fc=0 fc=0 

Whereas there are some results concerning the coefficients Bk, almost nothing 
is known about the coefficients Hk- Since the Zaremba problem is local, or 
better to say 'pseudo-local', all these coefficients have the form 



B2. = Jb^^^ + Jb^r + JbT + Jb'^Z (80) 

M El S2 So 

B2k^i = / 4^+1 + / 4^+1 + /d' (81) 

Si S2 So 



Hk = hk. (82) 
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Here the new feature is the appearance of the integrals over Eq, which com- 
phcates the problem even more, since the coefficients now depend on the 
geometry of the imbedding of the codimension 2 submanifold Sq in M that 
could be pretty complicated, even if smooth. 

Let us stress here that we are not going to provide a rigorous construction 
of the parametrix with all the estimates, which, for a singular boundary- 
value problem, is a task that would require a separate paper. For a complete 
and mathematically rigorous treatment the reader is referred to [^ and 
references therein. Here we keep instead to a pragmatic approach and will 
describe the construction of the parametrix that can be used to calculate 
explicitly the heat kernel coefficients B^ as well as Hk- 

First of all, we need to properly describe the geometry of the problem. 
Let us fix two small positive numbers 61,62 > 0. We split the whole manifold 
in a disjoint union of four different parts: 

= M^°* U M^""^ U M2^°^ U Mq^'^'^ . 

Here M^^*^ is defined as the set of points in the narrow strip M^^*^ of the 
manifold M near the boundary dM of the width 61 that are at the same 
time in a narrow strip of the width 62 near Sq 

M^""^ = {x G M I dist(x, dM) < 61, dist(x, Sq) < 62} . 

Further, M^'^'^ is the part of the thin strip M^^'^ of the manifold M (of the 
width 61) near the boundary dM that is near Si but at a finite distance from 
Eo, i.e. 

M^""^ = {x e M \ dist(x. El) < 61, dist(x, Sq) > 62} . 

Similarly, 

M^""^ = {x e M \ dist(x, S2) < 61, dist(x. So) > 62} . 

Finally, M™* is the interior of the manifold M without a thin strip at the 
boundary dM, i.e. 

M'°* = M\ [M]"""^ U M^""^ U M^""^) 
= {x e M \ dist(a;, dM) > 61} . 
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We will construct the parametrix on M by using different approximations 
in different domains. Strictly speaking, to glue them together in a smooth 
way one should use 'smooth characteristic functions' of different domains 
(partition of unity) and carry out all necessary estimates. What one has 
to control is the order of the remainder terms in the limit t — > and their 
dependence on ei and 62- Since our task here is not to prove the form of the 
asymptotic expansion (|79D , which is known, but rather to compute explicitly 
the coefficients of the asymptotic expansion, we will not worry about such 
subtle details. We will compute the asymptotic expansion as t — *> in each 
domain and then take the limit 81,82 — > 0. 

We will use different local coordinates in different domains. In M™* we 
do not fix the local coordinates; our treatment will be manifestly covariant. 

In M}'"'^ we choose the local coordinates as follows. Let 

{ei}, 

(2 = 2,... m), be the local frame for the tangent bundle TEi and x = (x*) = 
(x^, . . . , x™), (i = 2, . . . , m), be the local coordinates on Si. Let 

r = dist(x. Si) 

be the normal distance to Si ( 

r = 

being the defining equation of Si), and 

be the inward pointing unit normal to Si. Then by using the geodesic flow 
we get the local frame 

{N,e,} 

for the tangent bundle TM and the local coordinates x = (r, x) on M^^'^. 
The coordinate r ranges from to 81, 

0<r <8i. 

The local coordinates in M^'^'^ are chosen similarly. 

Finally, in M^^'^ we choose the local coordinates as follows. Let 
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(a = 3, . . . ,m), be a local frame for the tangent bundle TSq and let x = 
(£") = {x^, . . . , X™) be the local coordinates on Sq. To avoid misunderstand- 
ing we should stress here that now we use the same notation x to denote 
coordinates on Eq (not on the whole of dM). Let dist aM(a^, Sq) be the 
distance from a point x on dM to So along the boundary dM. Then define 

y = +distaM(a;, So) > 

if a; G Si and 

y = -distaM(a:, So) < 

if a; G S2. In other words, 

y = o 

on So ( 

r = y = 

being the defining equations of So), 1/ > on Si and 1/ < on S2. Let 

be the unit normal to Sq pointing inside Si. Then by using the tangential 
geodesic flow along the boundary (that is normal to Sq) we first get the local 
orthonormal frame 

{n{y,x),ea{y,x)} 

for the tangent bundle TdM. Further, let the unit normal vector field to the 
boundary N{y, x) be given. Then by using the normal geodesic flow to the 
boundary we get the local frame 

{N{r, y, x), n(r, y, x), ea(r, y, x)} 

for the tangent bundle TM and local coordinates {r,y,x) on Mq°*^. The 
ranges of the coordinates r and y are: 

< r < £1 

and 

-£2<y < £2- 

Finally, we introduce the polar coordinates 

r = pcos6, y = psinO . 
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To cover the whole Mq'^^ the angle 9 ranges from — 7r/2 to 7r/2 and p ranges 
from to some £3 (depending on Si and £2), 

< p < £3- 

The construction of the interior parametrix goes along the same lines as 
for manifolds without boundary (see lectures 2 and 4). For a finite ei > 
the diagonal of the heat kernel has the same asymptotic expansion as for 
manifolds without boundary. Therefore, by integrating over the interior part 
M™* and taking the limit ei ^ we find that the local interior coefficients 
h\. are the same as for manifolds without boundary in the smooth case, i.e. 

af , given by (§HP-(§g). 

The Dirichlet parametrix U^'^'^''^'^\t) in M"^^^ and Neumann parametrix 
f/bnd,(2)j^^-j in M}^"^^ are constructed along the same lines as the parametrix of 
a smooth boundary- value problem described in lecture 5. For finite ei, £2 > 
the diagonal of the parametrix has the same kind of asymptotic expansion 
as t — > with coefficients being homogeneous functions of r / \/t 



<^'"w~E''^c/r-'"(] 

j,_n V 



Vt'^'v^' 



fc=0 

After integrating the diagonal over M]'^'^ and taking the limit £1, £2 ~^ the 
contribution of all regular terms will vanish. The singular terms, which be- 
have like distributions near Ej, will give the Sj contributions to the boundary 
heat kernel coefficients 6^ . As the result we get the coefficients b]. in the 
form 






n=0 



xtrv^f/ril(l|e,x;e,x). (83) 

These are the standard boundary heat kernel coefficients for smooth Dirichlet 
and Neumann boundary conditions. They are listed for example in P^ , 2^ 
up to fc = 4. The first two have the form 



,(1),1 _ Jl),2_(. 
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6«'i = -6«'^ = -(47r)-('"-^)/Mimvl, 

ftW'i = 6^'^'' = (47r)-™/2dimy-ir, (84) 

where K is the trace of the extrinsic curvature (second fundamental form) 
of the boundary. 

The most comphcated (and the most interesting) is the case of the mixed 
parametrix in M^^^ since here the basic problem with frozen coefficients 
on So is two-dimensional. More precisely, in M^'^'^ the basic problem is on 
the half-plane. Since the origin is a singular point, we will work in polar 
coordinates introduced above. One can still use the scaling device described 
above. Since now we are working in the vicinity of the submanifold Sq, the 
coordinate p should also be scaled, i.e. 



p -^ ^p, 

t -^ eH, 

and similarly for the coordinates x', r' and p' . Then one needs to expand 
in e and develop a perturbation theory, which gives a recursion system that 
determines all coefficients Bk and Hk- The order at which the logt terms 
show up depends on the dimension of the manifold. 

We will restrict ourselves to the zeroth order of this perturbation theory 
only. This is enough for computation of the coefficient Bi. So, we are going 
to solve only the basic problem for operators with frozen coefficients at a 
point Xo on Sq. We choose normal coordinates on Eq at this point (with 
QabiS^, G, xq) = Sab)- Then the zeroth order operator Fq has the form 

Fo = -dl - -d, - \dl - d' (85) 

P P 



and the zeroth order normal is 



P 
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61 



Nn 



S2 P 



(86) 



Now the boundary operator is discontinuous, and there is a singularity at 
the origin p = 0. 

By separating the 'free' semiclassical factor due to Eq we get the ansatz 



u^'''^''\t\p,e,x,p\e\x') = (47rt)-(™-2)/^ 



X exp 






.'1 2 



4t 



^it\p,9;r',9'), (87) 



where \l/(t|p, ^; p', ^') is a two-dimensional heat kernel determined by the 
equation 

(d, - dl - -^d, - ^d^^ ^{t\p, 6; p', e') = 0, 



the initial condition 



1 



^(0+\p,9;p',9') = -S{p-p')S{9-9'), 
P 



the boundary conditions 



^{t\p,9;p',9' 



0, 



de^it\p,9;p',9'] 



and the asymptotic condition at infinity 



0, 



lim ^(t|p,^;p',^') = 0. 

Clearly, the heat kernel is also symmetric 

^it\p,9;p',9') = ^it\p',9';p,9). 

This problem can be solved by separating variables, employing the Hankel 
transform in the radial coordinate and evaluating a certain spectral series 
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of Bessel functions. As a result, we obtain the mixed leading parametrix in 
U^'''''^'\t\p,9;p',9') = Lit\p,9,x;p',9',x') 



M"^^^ 



+L(t|p, 9, x; p', -9' - TT, x') 
where 
L{t\p,9,x;p',9',x') = (47rt)-'"/2 

\x — x'p + p^ + p'^ — 2pp' cos( 



X exp 



X erf < \ — cos 



4t 



. t V 2 

with 

z 

erf iz) = —= I due "^ 
V TT J 



being the error function. 

The diagonal of the mixed parametrix is easily found to be 



^dtjrW = (4vrt)-™/^ 1 



• (Q\ ( P^COS^0\ ( p 

— Sign (6') exp I I — errc I —jn 

( p^cos^9\ , f p , ,, 
+ sign [9) exp I j errc I —= \sm9\ 

where sign (x) is the sign function, i.e. sign (x) = 1 for a; > and sign {x) = 
— 1 for X < 0, and 

erfc (z) = 1 — erf (z) 

POO 

= 27r-^/2 / due-''' (89) 
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is the complementary error function. 

Finally, we compute the integral of the diagonal of the mixed parametrix 
over Mq"^*^ and take the limit £1,62 —^ 0. We have 

Jm^ / tr yUrS\t)=}l^ J I dp P I d9tryUrS\t)^ 

Aft'"d So -7r/2 



for some finite £3 > ^/ef + e\ > 0. By computing the integrals and taking 
the limits we obtain 

^^lim^ f trvf/dtg;[,°^(t) = -t(2-"^)/2(47r)^"^/2^dim\/vol(So). 



jVjbnd 



l(2) 



This gives exactly the coefficient 63 i^i the heat trace asymptotic expansion 
(H), i.e. 

6f = -(47r)-(™-2)/2dimr — . (90) 
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Note Added in Proof 

The Zaremba problem considered in the section 6 was studied recently by 
Seeley in [^| . It has been shown that the logarithmic terms in the expansion 



([79| ) are absent, i.e. H^ = for any k, which confirms the conjecture of 



5B| . Seeley has also shown that the correct setting of the Zaremba problem 
involves an additional boundary condition along the singular set Sq. The 
general solution as well as the heat kernel asymptotics do depend on this 
additional condition. Our solution corresponds to the choice of most regular 
eigenfunctions close to Sq. Other solutions contain integrable singularities 
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